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Abstract 

We report on a search for the recently reported X(3872) state using 15.1 fb _1 of e + e~ data taken 
in the y/s = 9.46-11.30 GeV region. Separate searches for the production of the X(3872) in untagged 
77 fusion and e + e~ annihilation following initial state radiation are made by taking advantage of 
the unique angular correlation between the leptons from the decay J/ip l + U in X(3872) decay to 
tt + tt~ J /if;. No signals are observed in either case, and 90% confidence upper limits are established 
as (2J + l)r 77 (X(3872))^(X -► n+ir- J/Y>) < 12.9 eV and T ee (X(3872))B(X -► vr+vr-J/^) < 8.3 
eV. 



PACS numbers: 12.39.Mk,13.66.Lm,13.25.Gv,14.40.Gx 
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The Belle Collaboration recently reported the observation of a narrow state, X(3872), 
in the decay B ± — >K ± X, X — > 7i + 7i~J/ip, J/ijj — > l + l~ {I = e,fi) 0]. The observation was 
confirmed by the CDF II 0], D0 0], and BaBar 0] collaborations, with consistent results, 
M{X) = 3872 ± 1 MeV/c 2 , and T{X) < 3 MeV/c 2 . 

Many different theoretical interpretations of the nature of the X(3872 ) st ate and its 
possible quantum numbers have been proposed JEM 0, H, EH EH 2 IS 3 3- These 
include that (a) X(3872) is a charmonium state p,|6|; (b) X(3872) is a D°D*° loosely bound 
"molecular" state 00 since its mass is close to (M D o + M s *o) = 3871.3 ± 1.0 MeV/c 2 
(lif : and (c) X(3872) is an exotic state 0]. 

No positive signals for X(3872) have been observed in searches for the decay channels 
X(3872) -> 7Xci Q, 7Xc2, 7^M ttVJ/V 1>|J/^ Q, D+D~, D°D°, and D^tt 0, 
or for possible charged partners of X(3872) j20|. Yuan, Mo, and Wang [21] have used 22.3 
pb" 1 of BES data at ^/s = 4.03 GeV to determine the upper limit of r ee (X(3872))£(X -> 
ii + 7i~J/ib) < 10 eV (90% C.L.) for X(3872) production via initial state radiation (ISR). 
Belle pJZj has recently reported a small enhancement in the 7r + n~n° J/ip effective mass near 
the X(3872) mass. 

The variety of possibilities for the structure of X(3872) suggests that, irrespective of the 
models, it is useful to limit the J PC of X(3872) as much as possible. The present investigation 
is designated to provide experimental constraints for the J PC of X(3872) by studying its 
production in 77 fusion and ISR, and its decay into n + n^J/tp. Production of X(3872) in 
77 fusion can shed light on the positive charge parity candidate states, charmonium states 
2 3 P , 2 3 P 2 and 1 1 D 2 |^, and the 0~ + molecular state 0,0. ISR production can address 
the 1 vector state. 

The data used for this X(3872) search were collected at the Cornell Electron Storage 
Ring (CESR) with the detector in the CLEO III configuration The detector is cylindri- 
cally symmetric and provides 93% coverage of solid angle for charged and neutral particle 
identification. The detector components important for this analysis are the drift chamber 
(DR), Csl crystal calorimeter (CC), and muon identification system (MIS). The DR and 
CC are operated within a 1.5 T magnetic field produced by a superconducting solenoid 
located directly outside of the CC. The DR detects charged particles and measures their 
momenta and ionization energy loss (dE/dx). The CC allows precision measurements of 
electromagnetic shower energy and position. The MIS consists of proportional chambers 
placed between layers of the magnetic field return iron to detect charged particles which 
penetrate a minimum of three nuclear interaction lengths. 

The data consist of a 15.1 fb -1 sample of e + e~ collisions at or near the energies of the 
T(nS) resonances (n = 1-5), and in the vicinity of the AbA& threshold. Table [I] lists the six 
different initial center-of-mass energies and integrated luminosities of the data samples. 

Resonance production by untagged 77 fusion and by ISR have similar characteristics. 
The undetected electrons in untagged 77 fusion and the undetected radiated photons in 
ISR have angular distributions sharply peaked along the beam axis. Both processes have 
total observed energy (E tat ) much smaller than the center-of-mass energy, yfs, of the original 
e + e~ system, and have small observed transverse momentum. The detailed characteristics 
for 77 fusion and ISR resonance production are studied by generating signal Monte Carlo 
(MC) samples using GEANT 3.21/11 [23] to simulate the CLEO III detector. For X(3872) 
production by 77 fusion the formalism of Budnev et al. 24|| is used. For ISR resonance 



production the formalism of M. Benayoun et al. [25] is used. 

A fully reconstructed event has four charged particles and zero net charge. All 
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TABLE I: Data samples and MC determined detection efficiencies used for the present A(3872) 
search, {^/sl ) are the average center-of-mass energies of T(15 — 5S) and AjAf, threshold measure- 
ment and Ci(e + e~) is the e + e~ integrated luminosity at yfsi. The efficiencies e 77i j and eisR,i are 
the sums of the efficiencies t ee ^ and e W) i for electron and muon detection, respectively. The 77 
fusion/ISR separation, as described in the text, is applied to the respective MC samples. 





<v^> 


A(e + e~) 


77 Fusion 


ISR 




(GeV) 


(fb- 1 ) 












eiSR,i 


T(LS) 


9.458 


1.47 


0.128(4) 


0.160(4) 


0.288(6) 


0.065(3) 


0.083(3) 


0.148(4) 


T(2S) 


10.018 


1.84 


0.121(3) 


0.151(4) 


0.272(5) 


0.054(2) 


0.062(3) 


0.116(4) 


T(3S) 


10.356 


1.67 


0.115(3) 


0.137(4) 


0.252(5) 


0.042(2) 


0.043(2) 


0.085(4) 


T(45) 


10.566 


8.97 


0.123(4) 


0.145(4) 


0.268(6) 


0.0186(14) 


0.0165(13) 


0.0351(19) 


T(55) 


10.868 


0.43 


0.113(3) 


0.139(4) 


0.252(5) 


0.0025(5) 





0.0025(5) 


A^Aft threshold 


11.296 


0.72 


0.104(3) 


0.126(4) 


0.230(5) 


0.0001(1) 





0.0001(1) 



charged particles must lie within the drift chamber volume and satisfy standard require- 
ments for track quality and distance of closest approach to the interaction point. Events 
must also have detected E tot < 6 GeV. The X(3872) resonance corresponds to AM = 
M(7r + 7r-/ + /-) - M(Z+r) = 0.775 GeV/c 2 , and we designate AM = 0.63-0.7 and 0.85-0.92 
GeV/c 2 as background regions. Signal-to-background studies are performed to optimize 
signal efficiency and background suppression. Selection criteria optimized the efficiency for 
reconstructing 77 fusion MC events. The selection variables optimized are the total neutral 
energy (E neu ) of the event, total transverse momentum of the four charged tracks (pt r ), 
lepton pair invariant mass (M(/ + / - )) of the J/ip — > l + l~ decay, and particle identification 
(PID) of the charged tracks. Based on the optimization studies, events are selected with 
Ene« < 0.4 GeV and p tr < 0.3 GeV/c. Events with a J/^> -> e + e~ decay require both 
electron candidates to satisfy dE/dx and shower energy criteria consistent with the electron 
hypothesis, and to have invariant mass in the range M(e + e~) = 2.96-3.125 GeV/c 2 . Events 
with a J/ip — > decay require both muon candidates to appear as minimum ionizing 

particles in the CC, with at least one muon penetrating the number of interaction lengths in 
the MIS consistent with its momentum, and to have invariant mass in the range M(/i + /i~) 
= 3.05-3.125 GeV/c 2 . Each of the two pions recoiling against the J/ip is required to satisfy 
the dE/dx pion hypothesis. 

Figure ^ shows the AM distribution for data events which pass the selection criteria and 
have AM = 0.514-0.850 GeV/c 2 . A ip(2S) signal is clearly visible while no enhancement 
is apparent in the X(3872) region. The observed number of ip(2S) events is determined by 
fitting the ip{2S) region with a mass-independent background and a resonance whose shape 
is determined by fitting the ip(2S) peak in the ISR MC simulation. The observed number 
of ip(2S) is N JSR (ip(2S)) = 206 ± 15 events. A MC simulation predicts N ISR (ip(2S)) = 226 
±11 events. 

At \/s ~ 10 GeV, a feature unique to the ISR mediated production of a vector resonance 
which decays via 7T + 7r~J/il), J/ip — ► Z + /~ is the correlation between the angles 9i+ and 
61- in the laboratory system. Figure 121 shows the MC prediction for the two-dimensional 
cos(#) distributions for leptons from V(3872) decay for the ISR mediated and 77 fusion 
productions. As shown in Figure |2l a parabolic cut applied to the two-dimensional cos(#) 
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FIG. 1: Data events as function of AM = M(7r + 7r"/+/") - M(l + l~). The ip(2S) is clearly visible 
and no apparent enhancement is seen in the X(3872) region. 



3850804-002 




FIG. 2: MC predictions for the two-dimensional lepton pair cos(#) distributions for the X(3872): 
ISR (left) and 77 fusion (right). The lines indicate how the ISR resonance and 77 fusion samples 
are separated. 



distribution efficiently separates the events from the two production processes. With this 
cut, the 77 region contains a 0.6% contamination from ISR production, and the ISR sample 
contains a 14% contamination from 77 fusion production if we assume for an illustrative 
purpose that (2 J + 1)T 77 (X) = T ee (X). Here J is the total spin and T 77 (T ee ) is the 
two-photon (e + e~) partial width of X(3872). 

The efficiencies as determined by MC simulations of X(3872) production and decay fol- 
lowing 77 fusion and ISR are listed in Table HI The X(3872) and J/ip are decayed according 
to phase space in the MC simulations. The same selection criteria are applied to both MC 
samples except for the lepton pair cos(#) cut described above. 

The separate AM distributions for the data in the X(3872) search region for 77 fusion 
and ISR mediated resonance production are shown in Figure El The number of observed 
X(3872) events (iV 77j / > s7j(X(3872))) is determined by maximum likelihood fits of the AM 
data using mass-independent backgrounds and the appropriate detector resolution functions 
for the two production processes. The detector resolution functions are determined by the 
MC simulations fitted with double Gaussians which are illustrated in Figure H3 The 90% 
confidence upper limits on the observed number of X(3872) events in untagged 77 fusion 
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FIG. 3: Distributions of data events as function of AM = M(7r + 7r~Z + Z~) — M(Z + Z~) for 77 fusion 
(top) and ISR (bottom) events in the region AM = 0.7-0.85 GeV/c 2 . The mass resolution functions 
determined from MC simulations are shown on an arbitrary scale at AM = 0.775 GeV/c 2 . 



and ISR mediated resonance production are determined to be X 77i / S t ? (X(3872)) < 2.36 for 
both processes. 

The cross section for 77 fusion or ISR mediated production of the X(3872) resonance 
with total angular momentum J, and decay through 7i + 7i^J/ip, J/if> — > is 



' M 



)(2J+ l)T 11)ee (X)B{X -> n + n-J/ij) = 
X 77i/5fi (X(3872)) 



B(J/ip^l+l )EiA(e+e )e 11 j S R^{^/sl) 11 ,isR 

where C^jsr are constants, M = 3872 MeV/c 2 , ^/sT, £j(e + e~), and e^jsRj are as listed in 
Table EJ and &(y/si) T y,isR are as shown in Figure |U The branching fraction B(J/if> — ► l + l~) 
=_[5.91±0.07)% is the average PDG branching fraction of J /if) — > e + e~ and J /if) — > 
[16]. This leads to the 90% confidence upper limits 

(2 J + l)r 77 (X(3872))£(X -> Tr+Tr" J /if}) < 10.9 eV 

for X(3872) production in 77 fusion, and 

r ee (X(3872))£(X -> 7T+71- J/ifj) < 7.3 eV 

for X(3872) production via ISR. 

Systematic uncertainty in the above limits arises from possible biases in the detection 
efficiency and estimated background level. These are studied by varying the track quality, 77 
fusion/ISR separation, and selection criterion optimized in the signal-to-background studies. 
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FIG. 4: Cross sections for e + e~ collisions at \fs to produce reduced CM energy, vV, for 77 fusion 
with V? = 3872 MeV and ISR with \/7' = 3872 MeV and V? = 3686 MeV. 

Other systematic uncertainties are from the e + e~ luminosity measurement and J/if> — > Z + /~ 
branching fractions. Adding these in quadrature, the total systematic uncertainties in 77 
fusion and ISR are 18.5% and 13.2%, respectively. A conservative way to incorporate these 
systematic uncertainties is to increase the measured upper limits by these amounts. This 
leads to the 90% confidence upper limits 

(2J + l)r 77 (X(3872))i3(X -> Tr+Tr" J/tjj) < 12.9 eV 

for X(3872) which has positive C parity, and 

r ee (X(3872))£(X -> 7T+7T- J/ip) < 8.3 eV 

for X(3872) being a vector meson with J PC = 1 . 

If B(B ± ^K ± X (3872)) w ^(^±^^±^(2^)) = (6.8±0.4)xl0^ 4 is assumed, we ob- 
tain 8{X -> 7i + 7i-J/ip) w 0.02 from both the Belle Q and BaBar Jresults. This leads to 
90% confidence upper limits 

(2J + l)r 77 (X(38 72)) < 0.65 keV, 

and 

r ee (X(3872)) < 0.42 keV. 

The (2 J+l)r 77 (X(3872)) upper limit is almost 1/4 the corresponding values for Xco and 
Xc2i but it is nearly 6 times larger than the prediction for the 1 1 D 2 state of charmonium [26]. 
The upper limit for r ee (X(3872)) is comparable to the measured electron width of ^(3770) 
and is about 1/2 that of ^(4040). We also note that the ratio N ISR (X(3872))/N IS R(i>(2S)) 
< 0.01 (90% C.L.). 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 
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